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The role of torsion in f(R) gravity is considered in the framework of metric-affine formalism. We 
discuss the field equations in empty space and in presence of perfect fluid matter taking into account 
the analogy with the Palatini formalism. As a result, the extra curvature and torsion degrees of 
freedom can be dealt as an effective scalar field of fully geometric origin. From a cosmological point 
of view, such a geometric description could account for the whole Dark Side of the Universe. 
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I. INTRODUCTION 



In the last thirty years, some shortcomings came out in the Einstein General Relativity (GR) and several investi- 
gations started in order to study if alternative approaches to gravitational interaction are possible and self-consistent. 
00 ' Such issues come from Cosmology and Quantum Field Theory. In the first case, the presence of Big Bang singularity, 
CNJ ' flatness and horizon problems [l[ led to the result that Standard Cosmological Model || , based on GR and Standard 
Model of particle physics, is inadequate to describe the Universe at extreme regimes. On the other hand, GR does 
' not work for a quantum description of spacetime. Due to this facts and to the lack of a Quantum Gravity theory, 
Q H , alternative theories of gravity have been pursued in order to attempt, at least, a semi-classical scheme where GR and 
j ■ its positive results could be recovered. 

A fruitful approach has been that of Extended Theories of Gravity (ETG) which have become a sort of paradigm 
l—— '■ in the study of gravitational interaction. They are essentially based on corrections and enlargements of the Einstein 
| theory. The paradigm consists in adding higher-order curvature invariants and non-minimally coupled scalar fields 
into dynamics resulting from the effective action of Quantum Gravity [3|, |J] . 

Other motivations to modify GR come from the issue to recover Mach's principle Q. This principle states that 
C*") , the local inertial frame is determined by some average of the motion of distant astronomical objects @, so that 
gravitational coupling can be scale-dependent and related to some scalar field. This viewpoint leads to assume a 
varying gravitational coupling. As a consequence, the concepts of "inertia" and equivalence principle have to be 
00 ■ revised 0,0,11,0. 

Furthermore, every unification scheme as Supcrstrings, Supcrgravity or Grand Unified Theories, takes into account 
effective actions where non-minimal couplings to the geometry or higher-order terms in the curvature invariants come 
out. Such contributions are due to one-loop or higher-loop corrections in the high-curvature regimes. In particular, 
^ . this scheme has been adopted in order to deal with quantization in curved spacetimes and, as a result, the interactions 
between quantum scalar fields, the gravitational self- interactions and the background geometry yield correction terms 
m | in the Hilbert-Einstein Lagrangian |3~0f | - 

Moreover, it has been realized that such terms are inescapable if we want to obtain the effective action of Quantum 
Gravity on scales closed to the Planck length Higher-order terms in the curvature invariants, as R 2 , R^Rij, 

R^ kh Rijf.fr, ROR, RO k R or non-minimally coupled terms between scalar fields and geometry, as <fi 2 R, have to be 
added to the effective Lagrangian of gravitational field when quantum corrections are considered. For example, one has 
to stress that such terms occur in the effective Lagrangian of strings or in Kaluza-Klein theories, when the mechanism 
of dimensional reduction is used [HI]. 

From a conceptual point of view, there would be no a priori reason to restrict the gravitational Lagrangian to a 
linear function of the Ricci scalar R, minimally coupled with matter (T^]. The idea that there are no "exact" laws 
of physics but that the Lagrangians of physical interactions are "stochastic" functions - with the property that local 
gauge invariances (i.e. conservation laws) are well approximated in the low energy limit and that physical constants 
can vary - has been taken into serious consideration - see Ref. [14| and references therein. 

Besides fundamental physics motivations, all these theories have acquired a huge interest in cosmology due to the 
fact that they "naturally" exhibit inflationary behaviors able to overcome the shortcomings of Standard Cosmolog- 
ical Model (based on GR). The related cosmological models seem very realistic and capable of matching with the 
observations fla. [la. flTll . 



Furthermore, it is possible to show that, via conformal transformations, the higher-order and non-minimally coupled 
terms always correspond to Einstein gravity plus one or more than one minimally coupled scalar fields [H, [ljj [2(| 
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l2ll . |22| . More precisely, higher- order terms always appear as a contribution of order two in the equations of motion. 
For example, a term like R 2 gives fourth order equations (23[, R OR gives sixth order equations dHHl]) R0 2 R gives 
eighth order equations [25[ and so on. By a conformal transformation, any 2nd-order of derivation corresponds to a 
scalar field: for example, fourth-order gravity gives Einstein plus one scalar field, sixth order gravity gives Einstein 
plus two scalar fields and so on [2lL [2t| . This feature results very interesting if we want to obtain multiple inflationary 
events since a former early stage could select "very" large-scale structures (clusters of galaxies today), while a latter 
stage could select "small" large-scale structures (galaxies today) [24| ■ The philosophy is that each inflationary era is 
connected with the dynamics of a scalar field. Furthermore, these extended schemes naturally could solve the problem 
of "graceful exit" bypassing the shortcomings of former inflationary models H3] ■ 

Recently, ETG are going also to play an interesting role to describe the today observed Universe. In fact, the 
amount of good quality data of last decade has made it possible to shed new light on the effective picture of the 
Universe. Type la Supernovae (SNela) [28|, anisotropies in the cosmic microwave background radiation (CMBR) 29], 
and matter power spectrum inferred from large galaxy surveys [30| represent the strongest evidences for a radical 
revision of the Cosmological Standard Model also at recent epochs. In particular, the concordance KCDM model 
predicts that baryons contribute only for ~ 4% of the total matter - energy budget, while the exotic cold dark matter 
(CDM) represents the bulk of the matter content (~ 25%) and the cosmological constant A plays the role of the so 
called "dark energy" (~ 70%) [31j. Although being the best fit to a wide range of data [32|], the ACDM model is 
severely affected by strong theoretical shortcomings [33| that have motivated the search for alternative models [HI, [35| . 

Dark energy models mainly rely on the implicit assumption that Einstein's GR is the correct theory of gravity 
indeed. Nevertheless, its validity at the larger astrophysical and cosmological scales has never been tested |36| . 
and it is therefore conceivable that both cosmic speed up and dark matter represent signals of a breakdown in our 
understanding of gravitation law so that one should consider the possibility that the Hilbert - Einstein Lagrangian, 
linear in the Ricci scalar R, should be generalized. 

Following this line of thinking, the choice of a generic function f(R) can be derived by matching the data and by 
the " economic" requirement that no exotic ing redients have to be added. This is the underlying philosophy of what is 
referred to as f(R) gravity [37l.[38l[39l.[40l l4ltl42l. \I§L lij. liR IZej l47j . In this context, the same cosmological constant 
could be removed as an ingredient of the cosmic pie being nothing else but a particular eigenvalue of a general class 
of theories 

However f(R) gravity can be encompassed in the ETGs being a "minimal" extension of GR where (analytical) 
functions of Ricci scalar are taken into account. 

Although higher order gravity theories have received much attention in cosmology, since they are naturally able to 
give rise to accelerating expansions (both in the late and in the early Universe) and systematic studies of the phase 
space of solutions are in progress [H, HH, [H, [53[, it is possible to demonstrate that f(R) theories can also play 
a major role at astrophysical scales. In fact, modifying the gravity Lagrangian can affect the gravitational potential 
in the low energy limit. Provided that the modified potential reduces to the Newtonian one on the Solar System 
scale, this implication could represent an intriguing opportunity rather than a shortcoming for f(R) theories (see for 
example H M, HE IH H| ) . 

Furthermore, a corrected gravitational potential could offer the possibility to fit galaxy rotation curves without the 
need of dark matter [59|, HE EH, H2] • In addition, it is possible to work out a formal analogy between the corrections to 
the Newtonian potential and the usually adopted dark matter models. In general, any relativistic theory of gravitation 
can yield corrections to the Newton potential [63j) which, in the post-Newtonian (PPN) formalism, could give rise to 
tests for the same theory [H, EH, EE [66| . 

In this paper, we want to face the problem to study /(i?) gravity considering also torsion. Torsion theories have 
been taken into account firstly by Cartan and then where introduced by Sciama and Kibble in order to deal with spin 
in General Relativity (see [671 ] for a review). Being the spin as fundamental as the mass of the particles, torsion was 
introduced in order to complete the following scheme: the mass (energy) as the source of curvature and the spin as 
the source of torsion. 

Up to some time ago, torsion did not seem to produce models with observable effects since phenomena implying 
spin and gravity were considered to be significant only in the very early Universe. After, it has been proven that spin 
is not the only source of torsion. As a matter of fact, torsion field can be decomposed in three irreducible tensors, with 
different properties. In 68j, a systematic classification of these different types of torsion and their possible sources is 
discussed. This means that a wide class of torsion models could be investigate independently of spin as their source. 

In principle, torsion could be constrained at every astrophysical scale and, as recently discussed, data coming from 
Gravity Probe B could contribute to this goal also at Solar System level |69j . 

In (70 . [7l| . a systematic discussion of metric-affine f(R) gravity has been pursued. In particular, the role of 
connection in presence of matter has been studied considering the various possible matter actions depending on 
connection. The main result of these papers has been the evidence that matter can tell to spacetime how to curve as 
well as how to twirl. 
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In this paper, following the same philosophy, we want to show that, starting from a generic f(R) theory, the 
curvature and the torsion can give rise to an effective curvature-torsion stress-energy tensor capable, in principle, to 
address the problem of the Dark Side of the Universe in a very general geometric scheme. We do not consider the 
possible microscopic distribution of spin but a general torsion vector field in f(R) gravity. 

The layout of the paper is the following. In Sees. II and III, we derive the metric-affine field equations of f(R) 
gravity with torsion in empty space and in presence of matter, respectively. Sec. IV is devoted to the discussion of the 
formal equivalence with scalar-tensor theories, while applications to Friedmann-Robertson- Walker (FRW) cosmology 
are discussed in Sec.V. Summary and conclusions are drawn in Sec. VI. 



II. FIELD EQUATIONS IN EMPTY SPACE 



Let us discuss the main features of a f(R) gravity considering the most general case in which torsion is present 
in a U4 manifold. In a metric-affine formulation, the metric g and the connection T can be, in general, considered 
independent fields. More precisely, the dynamical fields are pairs (g, T) consisting of a pseudo-Riemannian metric g 
and a metric compatible linear connection T on the space-time manifold M. The corresponding field equations are 
derived by varying separately with respect to the metric and the connection the action functional 



A(g,T) = J y/\g\f{R)da (1) 

where / is a real function, R(g, T) — g 1 ^ Rij (with R^ := R h ihj) is the scalar curvature associated with the connection 
r and ds := dx 1 A ■ • • A dx . Throughout the paper we use the index notation 

dT h r)V h 

nh _ jk _ UL ik , p hy, P _ y hy p /r>\ 

kij ~ dx' 1 dx j l P i k 3P ik ^ ' 

for the curvature tensor and 

d d 

V A^J = r */^T ( 3 ) 

for the connection coefficients. 

In order to evaluate the variation 8 A under arbitrary deformations of the connection, we recall that, given a metric 
tensor gij , every metric connection T may be expressed as 

T h = f h ~K h (A) 
y y y v > 

where (in the holonomic basis | — — , dx 1 1 ) T i4 h denote the coefficients of the Levi-Civita connection associated with 
v \dx l J 3 

the metric g^j, while -Ky' 1 indicate the components of a tensor satisfying the antisymmetry property K i j71 = —K i \ 

This last condition ensures the metric compatibility of the connection T. 

In view of this, we can identify the actual degrees of freedom of the theory with the (independent) components of 

the metric g and the tensor K. Moreover, it is easily seen that the curvature and the contracted curvature tensors 

associated with every connection Q can be expressed respectively as 



R iqj — R h iqj + ^ ' jK qi h — VqKji + K-^K j 1 — K - P K, h (5a) 



and 



% = R h ihj - Rij + VjKj 1 - V h K 3l h + K H p K hp h - KjK Jp h (5b) 



where R h i q j and Rij = R ! j k , are respectively the Riemann and the Ricci tensors of the Levi-Civita connection f 

associated with the given metric g, and V indicates the Levi-Civita covariant derivative. 

Making use of the identities (|5bj) , the action functional (p} can be written in the equivalent form 



A(g, T) = J V\g\f(g 13 (Rij + VjK hi h - V h K Si h + K/K hp h - KjK Jp h )) ds (6) 
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more suitable for variations in the connection. Taking the metric g fixed, we have the identifications STy = 5K i3 h 
and then the variation 



pats h 
tip 



5A = J ^\f{R)g^(V 3 5K hl h - V fc Mf/ + 6K/K hp h + K/SK, 

-5K h p K ]p h -K h p 5K Jp h )ds 



(7) 



Using the divergence theorem, taking the antisymmetry properties of K into account and renaming finally some 
indexes, we get the expression 



SA = 



+fK Jt h ] SK^ds 



The requirement 8 A = yields therefore a first set of field equations given by 

K p S h — K p S h - K h + K h = — — (S p S h -6 p S h ) 
pj 1 pi 3 v ""ji f'dxP ^ 3 

Considering that the torsion coefficients of the connection T are T.-,, h :— r,-,, /l — T^ h = —K 
to antisymmetry) T p p = —K pi p , eqs. ([9]) can be rewritten as 



or, equivalently, as 



h = — L^Zl (s p S h -6 p S h 



(8) 
(9) 

K ji h and thus (due 
(10) 

(11) 



In order to study the variation SA under arbitrary deformations of the metric, it is convenient to resort to the 
representation (TTJ). Indeed, from the latter, we have directly 



SA 



thus getting the second set of field equations 



f'(R)R ij - ^f(R)9ij 



j ij ds 



f'(R)R m - ^f(R)g t3 = 



(12) 



(13) 



Of course, one can obtain the same equations (fT3f starting from the representation ((6]) instead of {T]). In that case, 
the calculations are just longer. 

As a remark concerning eqs. (|13l) . it is worth noticing that any connection satisfying eqs. ([4]) and (fTTj) gives rise 
to a contracted curvature tensor automatically symmetric. Indeed, since the tensor K coincides necessarily with 
the contorsion tensor, namely 



from eqs. ([TTj) we have 



being 



K - = - l-T 



3 i 1 13 1 



K l0 h = \{T^-T p g ph g t3 ) 



Ti := T ih h 



js_df_ 

'if dx i 

Inserting eq. (| 15(1 in eq (|5bjl . the contracted curvature tensor can be represented as 

Rij = Rij + ^jTi + \v h T h g l0 + ^T t T 3 - ^T h T h 9ij 



(14) 



(15) 



(16) 



(17) 
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The last expression, together with eqs. (fT6|) , entails the symmetry of the indexes i and j. Therefore, in eq. (fl3|) we 
can omit the symmetrization symbol and write 

f'(R)Rij - \f{R)9n = o (18) 



Now, considering the trace of the equation (fT8|) . we get 

f'(R)R-2f(R) = (19) 

The latter is an identity automatically satisfied by all possible values of R only in the special case f(R) = aR 2 . In 
all other cases, eq. ([T9| represents a constraint on the scalar curvature R. 

As a conclusion follows that, if f(R) ^ aR 2 , the scalar curvature R has to be constant (at least on connected 
domains) and coincides with a given solution value of (|19[) . In such a circumstance, eqs. (lilt imply that the torsion 
T 4j h has to be zero and the theory reduces to a /(i?)-theory without torsion. In particular, we notice that in the 
case f(R) — R, eq. (TT9"]) yields R = and therefore eqs. (fT%|) are equivalent to Einstein's equations in empty space 
Rij — 0. On the other hand, if we assume f(R) = aR 2 , we can have non-vanishing torsion. In this case, by replacing 
eq. p^|) in eqs. (fTTj) and (Tig]) , we obtain field equations of the form 

Rij - -Rga = (20a) 

T . h -- ——d h + ——S h (20b) 

Finally, making use of eq. (fT7|) and the consequent relation 

R = R + 2V h T h - -T h T h (21) 
3 

in eqs. (|20p . we can separately point out the contribution due to the metric and that due to the torsion. In fact, 
directly from eqs. (|20a|) we have 

Rij - \Rtoi = - l^jTi + \v h T h 9l] - 2 -T{T 3 + ±T h T h 9i j (22) 

while from the "trace" T l := T lh h = -^§§ of eqs. ((20b)) . we derive 

^-{r + 2VhT h - ^T h T h ^j =-!(&+ 2V h T h - ^T h T h ^j T t (23) 

Eqs. (|2"2"|) and (|2^|) are the coupled field equations in vacuum for metric and torsion in the f{R) = aR 2 gravitational 
theory. 

III. FIELD EQUATIONS IN PRESENCE OF MATTER 

The presence of matter is embodied in the action functional (jTJ) by adding to the gravitational Lagrangian a suitable 
material Lagrangian density C rn , namely 

A(g,T) = J (y\g]f(R) + An) ds (24) 

Throughout the paper we shall consider material Lagrangian density C m not containing terms depending on torsion 
degrees of freedom as in [7lj |. The physical meaning of this assumption will be discussed later. In this case, the field 
equations take the form 

/'(«)«« - \nmsu = S.j (25a) 
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1 SH 

where Ey := ——j=-^j plays the role of the energy-momentum tensor. From the trace of eq. (|25a|) . we obtain a 

fundamental relation between the curvature scalar R and the trace £ := g lJ Y>ij, which is 

f'(R)R-2f(R) = E (26) 

(see also [72j and references therein). In what follows, we shall systematically suppose that the relation (f26l) is 
invertible and that £ ^ const., thus allowing to express the curvature scalar R as a suitable function of S, namely 

R = (27) 

With this assumption in mind, using eqs. (|2"6"|) and (|2"?| we can rewrite equations (|25a[) and (|25b[) in the form 

v-- w ^) a ™ w?-*» (28b » 

Moreover, making use of eqs. (|17|) and (l2Tj) . in cq. (|28a|) we can decompose the contracted curvature tensor and the 
curvature scalar in their Christoffcl and torsion dependent terms, thus getting an Einstein-like equation of the form 



+l^hT h 9ij ~ \TiTj - \T h T h 9lJ 



(29) 



Now, setting 



<P ■= f(F(Z)) (30) 



from the trace of eqs. (|28b[) . we obtain 



3 dip 



Ii:=T *--^ (31) 
Therefore, substituting in eqs. (|29|) . we end up with the final equations 

6 1 ~ 1 1 / 3 dip dip ~ dip 3 dip dip hk 

*Ui - 2 R 9ij ~ ~Ui + ^ [~2d^d^I + d^ + Adx^dx* 9 9ij 



(32) 



where we defined the effective potential 



V&) ~ \ {ipF-\{f)-\ip)) + <f?(f')-H<P)] (33) 

Eqs. (J22) may be difficult to solve, neverthless we can simplify this task finding solutions for a conformally related 
metric. Indeed, performing a conformal transformation of the kind g~ij = ipgij, eqs. (|32[) may be rewritten in the 
easier form (see, for example, [72l. 175ll76| ) 

1-11 

tki ~ 2 R 9ij = -Sij - -^ v VP)9ij (34) 

where Rij and R are respectively the Ricci tensor and the Ricci scalar curvature associated with the conformal metric 

9i3 ■ 

Concerning the connection T, solution of the variational problem 8 A = 0, from eqs. ([!]), (fT5j) and (pITj) . one gets the 
explicit expression 



r ,l -f h + -L^6 h -—^ a v h a ■ (35) 



7 



We can now compare our results with those obtained for f(R) theories in the Palatini formalism [44], [70, O, \T% 
IzE Iz3) IzH ■ If both the theories (with torsion and Palatini-like) are considered as "metric" , in the sense that the 

dynamical connection T is not coupled with matter f ^™ — O^j and it does not define parallel transport and covariant 

derivative in space-time, then the two approaches are completely equivalent. Indeed, in the "metric" framework, the 
true connection of space-time is the Levi-Civita one associated with the metric g and the role played by the dynamical 
connection T is just to generate the right Einstein-like equations of the theory. Now, surprisingly enough, our field 
equations (|32[) are identical to the Einstein-like equations derived within the Palatini formalism [751 ] . 

On the other hand, if the theories are genuinely metric-affine, then they are different even though the condition 

5 Cm 

— — = holds. In order to stress this point, we recall that in a metric-affine theory the role of dynamical connection 
ST 

is not only that of generating Einstein-like field equations but also defining parallel transport and covariant derivative 
in space-time. Therefore, different connections imply different space-time properties. This means that the geodesic 
structure and the causal structures could not obviously coincide. For a discussion on this point see [78| . Furthermore, 
it can be easily shown that the dynamical connection (|35[) differs from that derived within the Palatini formalism. 
Indeed the latter results to be the Levi-Civita connection f associated with the conformal metric g — <pg [74], [75| , 
while clearly (|35[) is not. More precisely, (|55"|) is related to T by the projective transformation 

' - =r « ' (36) 

which is not allowed in the present theory because, for a fixed metric <?, the connection (|36[) is no longer metric 
compatible. 

To conclude, we notice that eqs. (|34p are deducible from an Einstein-Hilbert like action functional only under 
restrictive conditions. More precisely, let us suppose that the material Lagrangian depends only on the components 
of the metric and not on its derivatives as well as that the trace S = Eyg y is independent of the metric and its 
derivatives. Then, from the identities 

/FT 2 /FT d 1 9 A v. 1 SLm 1 3C m 

V\9\=<P Vwh 7rTr = _ 7F77 and E ^ = TfrTTT = /TiTTF ( 37 ) 

og J V og J V\9\ ° g V Iff I " 9 

we have the following relation 

= ~^fU""7FF : = V^-v ( 38 ) 



dg l - 

In view of this, and being tp = y(E), it is easily seen that eqs. (|34[) may be derived by varying with respect to g 13 the 
action functional 



Mg) 



ds (39) 



Therefore, under the stated assumptions, f(R) gravity with torsion in the metric framework is conformally equivalent 
to an Einstein-Hilbert like theory. 



IV. EQUIVALENCE WITH SCALAR-TENSOR THEORIES 



The above considerations directly lead to study the relations between f(R) gravity with torsion and scalar-tensor 
theories with the aim to investigate their possible equivalence. To this end, we recall that the action functional of a 
(purely metric) scalar-tensor theory is 



1.91 



<pR ifiip 



(40) 



where ip is the scalar field which, depending on the sign of the kinetic term, could assume also the role of a phantom 
field [73|, ipi := — — and UUp) is the potential of tp. For UUp) = such a theory reduces to the standard Brans-Dickc 

theory [5j. The matter Lagrangian £ m (9iji "0) i s a function of the metric and some matter fields if>; uiq is the so called 
Brans-Dicke parameter. The field equations derived by varying with respect to the metric and the scalar field are 



His 2 H9ij - p ij + ^ 



1 h 

ViPj - ~<Ph<P 9ij 



- {^jVi - ^h<p h gij) - -^-9ij 



(41) 
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and 



2uJo V h( p h + R-^ hV h -U> = (42) 



p ip 

where Si, := == = m . and U 1 := -r— . 

Taking the trace of eq. (j4"Tj) and using it to replace R in eq. (|42|) . one obtains the equation 

{2uj + 3) V h ip h = S + <pU'- 2U (43) 

2 

By a direct comparison, it is immediately seen that for wo = — § and UUp) — —VUp) (where VUp) is defined as in 

p 

eq. ([53"| ). eqs. (|41 [) become formally identical to the Einstein-like equations (J32) for a /(i?) theory with torsion. 
Moreover, in such a circumstance, eq. (|4"3")l reduces to the algebraic equation 

E + tpU' - 2U = (44) 

relating the matter trace E to the scalar field </?, exactly as it happens for f(R) gravity. In particular, it is a 
straightforward matter to verify that (under the condition /" ^ 0) eq. (|44|) expresses exactly the inverse relation of 
(f3"0")) . namely 

S = F -i ((/r i M) ^ V=/'(F(E)) (45) 
being F-^X) = f(X)X - 2f(X). In fact we have 

U(<p) = lv(y,) = 1 [F^tf/TV)) + ^(/T 1 ^)] = ^(Z')- 1 ^) - W)-H<P))] (46) 



so that 

_^ <P _ f e ^-i 

/"((/')-%)) /"((/O-Hv)) 



and then 

S = -ifiU'ifp) + 2U(<p) = n(fT 1 i<PW)- 1 (<P) - 2/((/')- 1 M) - F^C/TV)) (48) 

As a conclusion follows that, in the "metric" interpretation, f(R) theories with torsion are equivalent to wo = — § 
Brans-Dicke theories. 

Of course, the above statement is not true if we regard f(R) theories as genuinely metric-affine ones. Nevertheless, 
also in this case it is possible to prove the equivalence between f(R) theories with torsion and a certain class of 
Brans-Dicke theories, namely luq — Brans-Dicke theories with torsion [761 ]. 

In this regard, let us consider the action functional 



A(g, r, <p) = J [y/\g\ (pR - U(ip)) + Cm\ ds (49) 

where the dynamical fields are respectively a metric gij, a metric connection and a scalar field p. As mentioned 
above, the action (|49p describes a Brans-Dicke theory with torsion and parameter ujq = 0. 
The variation with respect to ip yields the field equation 

R = U'(ip) (50) 

To evaluate the variations with respect to the metric and the connection we may repeat exactly the same arguments 
stated in the previous discussion for f(R) gravity. Omitting for brevity the straightforward details, the resulting field 
equations are 

V~5?!?(W-?*) <«u 

and 

R ij ~ l R 9ij = -Ey - ij-U(<p)gij (52) 
2 tp Zip 
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Inserting the content of eq. ([50)) in the trace of eq. (|52|) 

-E-~U(<p)+R=0 (53) 

ip ip 

we obtain again an algebraic relation between £ and ip identical to eq. (J3H). 

2 

Therefore, choosing as above the potential U(ip) = — VUp), from (|44|) we get ip = /'(F(E)). In view of this, 

f 

decomposing Rij and R in their Christoffel and torsion dependent terms, eqs. ([STj) and (|52|) become identical to 
eq. (f3"Tj) and (132[) respectively. As mentioned previously, this fact shows the equivalence between f(R) theories and 
luq = 0- Brans-Dicke theories with torsion, in the metric-affine framework. These considerations can be extremely 
useful in order to give a geometrical characterization to the Brans-Dicke scalar field. 



V. APPLICATIONS TO FRW COSMOLOGY 



We have seen that the field equations (|32|) may be recast in the form (|34f by performing a conformal transformation 
(jij — pgij. In order to apply the above considerations to FRW cosmological models, let us suppose that is the 
energy-momentum tensor of a cosmological perfect fluid with a negligible pressure and energy density p (dust case), 
namely 

E ij = pU l U j (54) 
where p = p(r) only depends on the cosmic time and U l is the four velocity of the fluid satisfying the condition 

QijWW = -1 (55) 

U' L 

From now on, we shall suppose ip > (a sufficient condition for this is /' > 0) so that the vector field U l := 

represents the four velocity of the fluid with respect to the conformal metric g^, while Ui := U J g~ji — J^pUi denotes 
the corresponding dual relation. In view of this, the identity 

— Sj,- = —UiUj = \%j (56) 
ip J p J ip 2 

holds, where we have defined = pUjJj. Consequently eqs. (|34[) may be rewritten as 

Ga = i hij - iv^Jfti) (57) 

where Gy is the Einstein tensor in the barred metric. We look for a FRW solution of {57]), being 

ds 2 = -dt 2 + a 2 {t) [dil> 2 + x 2 d9 2 + x 2 sin 2 6 dtp 2 ] (58) 

Therefore, once a solution g~ij is found, also the conformal metric gij = —g%j (solution of the starting equations (|32[> ) 
will be of the FRW form. Indeed, the line element associated to gij is 

ds 2 = L dt 2 + ^ to 2 + X 2 d9 2 + x 2 sin 2 9 d^ 2 } (59) 

¥>(*) V{t) 

so that, by performing the time variable transformation 

dr := i dt (60) 

it may be expressed as 

ds 2 = -dr 2 + A 2 {t) [di> 2 + x 2 dd 2 + x 2 sin 2 0dcf> 2 ] (61) 
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with A := . The field equations (l57|) reduce to the Friedmann equations 



and 



, / a\ 2 3k p V(tp) ,„„ , 



a f a\ 2 k V(p) ,„„, , 

2-+ - + — = -m (62b) 
a \a J or p a 



which are the cosmological equations arising from our theory. 

For the sake of completeness, let us derive the conservation laws of the theory. The Bianchi identities of eq. $5 
give 

J_E« _ \v{ V )gA = (63) 

where V denotes the covariant derivative with respect to the Levi-Civita connection associated with g~ij. In FRW 
metric, cqs. (|6"3")l reduce to the continuity equation 

which completes the cosmological dynamical system. 



VI. DISCUSSION AND CONCLUSIONS 



f(R) gravity seems a viable approach to solve some shortcomings coming from GR, in particular problems related 
to quantization on curved spacetime and cosmological issues related to early Universe (inflation) and late time dark 
components. Besides, the scheme of GR is fully preserved and f{R) can be considered a straightforward extension 
where the gravitational action has not to be necessarily linear in the Ricci scalar R. 

In this paper, we have discussed the possibility that also the torsion field could play an important role in the dynamics 
being the Ua manifolds the straightforward generalization of the pseudo-Riemannian manifolds V4 (torsionless) usually 
adopted in GR. 

As discussed above, torsion field, in the metric-affinc formalism, plays a fundamental role in clarifying the relations 
between the Palatini and the metric approaches: it gives further degrees of freedom which contribute, together 
with curvature degrees of freedom, to the dynamics. The aim is to achieve a self-consistent theory where unknown 
ingredients as dark energy and dark matter (up to now not detected at a fundamental level) could be completely 
" geometrized" . Torsion field assumes a relevant role in presence of standard matter since it allows to establish a 
definite equivalence between scalar-tensor theories and f(R) gravity, also in relation to conformal transformations. 

Finally, an important point deserves a further discussion in relation to the above results. Let us consider the 
cosmological equation (|62ap . In the Ihs, it is clear that standard matter p and the effective cosmological constant 

A-eff — pl a y two distinct role into the dynamics: their evolution is "tuned" by the scalar field ip (i.e. f'(R)). The 

first term could be relevant for large scale clustered structures (always involving baryonic matter and dark matter), 

the second term can be read as dark energy. If at present epoch they are ~ — ^— , this reveals a simple mechanism 

to explain why we are today observing Qm ~ 0.3 and Q\ ~ 0.7. Furthermore, if the field ip at denominator is small 
(that is f'{R) is small) this could be the reason why the amounts of dark energy and dark matter result huge today. 

As a toy model, let us take into account the well known f(R) = R + aR 2 theory where, obviously, f'(R) = l + 2<xR. 
As above, the matter stress-energy tensor is Ey = pUiUj and then Eq. ([26|) becomes 

(1 + 2aR)R - 2R - 2aR 2 = -p < — > R = p. (65) 

We have 



ip(p) = f'{R{p)) = l + 2ap 



(66) 
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and then the term becomes 

^ 2 



(1 + 2apf 

Let us consider now the potential term 



Being (/') 1 (ip) = p, one has 



and then we conclude that 



and 

V(<p(p)) ap 2 



ip 3 2(1 + 2ap) 



2 ' 



(67) 



V(<p) = \ [<pF- 1 ((f')- 1 {<P)) + V 2 (/')- 1 (¥')] • (68) 



l y{f)-\y) = 1(1 + 2apf±-{l + 2ap 1) = 1(1 + 2apf p , (69) 
and considering the relation F~ 1 {Y) = f'(Y)K - 2f(Y), it is 

l F -\{f)-\v)) = \F-\p) = -p. (70) 

We have also 

^- 1 ((/')- 1 M) = -^^, (71) 



VMp)) = + M (72) 



(73) 



These arguments show that the condition — ~ — — can be simply achieved leading to comparable values of Qm 
and Oa- A detailed discussion of these topic, also in relation with data, will be the argument of a forthcoming paper. 
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